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1. Introduction 

It is widely accepted that QCD undergoes a transition from the hadronic to the so-called 
"quark-gluon plasma" phase for large enough temperature, depending on the baryon chemical 
potential. Determining the shape and the nature of this transition (or pseudo-critical) line in the 
temperature - chemical potential plane is of central interest in cosmology, in astrophysics and in 
the phenomenology of heavy ion collisions. The lattice approach, which is the natural tool to 
face a non-perturbative problem like this, is plagued, however, by the well-known "sign problem": 
for non-zero chemical potential the QCD fermion determinant becomes complex and the standard 
Monte Carlo importance sampling is unfeasible. 

Several strategies have been invented to circumvent this problem (for a review, see ^ and [^). 
Here, we concentrate on one of these approaches, the method of analytic continuation, first used in 
Refs. 1^] and [Q]. The idea behind this method is very simple: numerical simulations are performed 
at imaginary chemical potential, }JL = ipLi, for which the fermion determinant is real, then Monte 
Carlo determinations are interpolated by a suitable function and finally this function is analytically 
continued to real values of }X. This method is rather powerful since temperature and chemical 
potential can be varied independently and there is no limitation from increasing lattice size, as in 
methods based on reweighting. There is, however, an important drawback: the periodicity of the 
partition function and the presence of non-analyticities arising for imaginary values of the chemical 
potential [||] limit the region useful for numerical determinations to the strip < Hj/T < n/3. This 
implies that the accuracy in the interpolation of the results at imaginary chemical potential has a 
strong impact on the extension of the domain of real jJ. values reachable after analytic continuation. 

Although the method is designed to infer the behavior of an observable with the real chemical 
potential from the knowledge of its dependence on the imaginary chemical potential, the idea 
to analytically continue the pseudo-critical line itself has been extensively applied [^, ^, ||] (see 
Ref. for a discussion on the reliability of this application of the method). 

A control on the systematics of the method of analytic continuation and possible insights for 
its improvement can be achieved by testing it in theories which do not suffer the sign problem, by 
direct comparison of the analytic continuation with Monte Carlo results obtained directly at real 



/I [Q, ^, |I0[ [TT|]. In Ref. [|TT|], in particular, a high-precision numerical analysis in SU(2) (or 2- 
color QCD) with = 8 degenerate staggered fermions has shown that, for temperatures above the 
pseudo-critical one at zero chemical potential, the extrapolation to real n improves considerably 
if ratio of polynomials are used instead of simple polynomials in interpolating the behavior with 
imaginary n of some test observables, this validating a proposal formulated in Ref. 



In this work we extend the numerical analysis of Ref. [ |11[ ] to the study of the analytic contin- 
uation of the pseudo-critical line. The strategy is the following: 

• for several fixed values of the chemical potential, both real and imaginary, we determine the 
(pseudo-)critical j8's by looking for peaks in the susceptibilities of a given observable; 

• we interpolate the determinations of the critical j8 's at imaginary chemical potential with an 
analytic function of /x, to be then extrapolated to real chemical potential; 

• we compare the extrapolation with the determinations of the critical jS's at real chemical 
potential. 
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Figure 1: (Left) Phase diagram in the {T, 9) plane according to Ref. (Right) Tentative phase diagram in 
the (r, 9) plane after the inclusion of the chiral pseudo-critical lines. 




Figure 2: Phase diagram in the (j3,/i/,/iff)-space; is the number of colors, A^^ the extension of the lattice 
in the temporal direction. 



The aim of this investigation is to verify the possibility to analytically continue the pseudo-critical 
line and, in affirmative case, to optimize the choice of the interpolating function to be used. For a 
better control of the systematics, we have repeated the outlined strategy for three different observ- 
ables (chiral condensate, Polyakov loop, plaquette). 

2. Theoretical background 

Long ago Roberge and Weiss have shown [^] that the partition function of any SU(A'^) gauge 
theory with non-zero temperature and imaginary chemical potential, /i = //X/, is periodic in 6 = 
Hi/T with period 2n/N and that the free energy F is a regular function of for T < Te, while it is 
discontinuous at = 2n{k + 1 /2)/N, k = 0,1,2, for T > Te, where Te is a characteristic tem- 
perature, depending on the theory. The resulting phase diagram in the {T, 0)-plane is given in Fig. ^ 



3 



Analytic continuation of the critical line in 2-color QCD at nonzero T and /i 



Alessandro Papa 




(left), where the vertical lines represent first order transition lines. This structure is compatible with 
the /I —pL symmetry, related with CP invariance, and with the Roberge-Weiss periodicity. The 
/I/ -dependence of any observable is completely determined if this observable is known in the strip 
< 6 < n/N. These predictions have been confirmed numerically in several cases, studying the 
behaviour of quantities like the Polyakov loop and the chiral condensate [|6l, ||, ^]. 

A phase diagram like that in Fig. [I] (left) would imply the absence of any transition along 
the T axis in the physical regime of zero chemical potential for any value of N, of rif and of the 
quark masses, which cannot be true. Therefore, it is necessary to admit that the phase diagram 
in the {T, 0)-plane is more complicated than in Fig. [I] (left). The simplest possibility is given in 
Fig. [l] (right), where the added lines generally represent transitions which can be first order, second 
order or crossover. The temperature Tc is the pseudo-critical one for the transition at zero chemical 
potential. In Fig. |2|the phase diagram of Fig. |l| (right) has been redrawn the (j3,/i/)-plane, where 
j8 = 2N/g^, pLi = ajJLi is the imaginary chemical potential in lattice units and it has been used the 
fact that r = 1/ {aNi), with A^^ the temporal extension of the lattice. In Fig. ^ also the /i«-axis has 
been included, with a sketch of the continuation of the pseudo-critical Une on the (j8,/i«) -plane. 
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Table 1: fi^ni as determined from the peak of the susceptibilities of chiral condensate, Polyakov loop and 
plaquette. 



A 


chiral condensate 


Polyakov loop 


plaquette 


0.30/ 


1.5097(33) 


1.5023(60) 


1.5048(62) 


0.20/ 


1.4548(42) 


1.4355(47) 


1.4476(94) 


0.15/ 


1.4258(51) 


1.441(23) 


1.418(10) 


0.10/ 


1.4228(55) 


1.4089(63) 


1.415(10) 


0. 


1.4072(35) 


1.394(17) 


1.4060(45) 


0.20 


1.3551(91) 


1.352(13) 


1.356(12) 


0.30 


1.271(16) 


1.267(26) 


1.286(15) 



Table 2: Parameters of the fit with a polynomial A + Bfi^ of the data at imaginary chemical potential; the last 
column gives the extrapolation of the pseudo-critical line at the imaginary chemical potential corresponding 
to the first RW transition line. 



observable 


A 


B 


;^^/d.o.f. 




chiral condensate 


1.4071(27) 


-1.140(50) 


0.91 


1.5829(81) 


Polyakov loop 


1.3931(55) 


-1.18(10) 


0.85 


1.575(17) 


plaquette 


1.4042(39) 


-1.104(83) 


0.48 


1.574(13) 



3. Numerical results 

We performed numerical simulations on a 16^ x 4 lattice of the SU(2) gauge theory with 
= 8 degenerate staggered fermions having mass am = 0.07, by means of a Hybrid Monte Carlo 
algorithm with dt = 0.01. Simulations have been performed on the APEmille crate in Bari and on 
the computer facilities at the INFN apeNEXT Computing Center in Rome. 

The observables we have considered are the chiral condensate, the Polyakov loop and the 
plaquette; for each of them we have looked for the peak in the susceptibility for varying j8, at 
some fixed values of the chemical potential, both real and imaginary. For each observable and for 
each fixed value of /t, we have taken a few data points for different /3 values (statistics ~ 20000) 
and smoothed out the susceptibility by the multi-histogram extension of the Ferrenberg-Swendsen 
reweighting method [pj[]. The uncertainty on the position of the peaks has been evaluated by the 
bootstrap method. 

In Table [l] we summarize our preliminary results for the critical j3 at each value of the chemical 
potential we considered. In Fig. ^ we show for illustration purposes the susceptibility of the chiral 
condensate at /t=0.15/, /t=0. and /i=0.30. 

We have then looked for an interpolation of the data of jScnt for imaginary and zero chemical 
potential (i.e. for the first 5 entries in Table and have repeated this procedure for each of the 
three observables considered. In all cases we have found that the optimal interpolating function is 
a polynomial of the form A +B/t^. If different functions are used, such as larger order polynomials 
or ratio of polynomials, the fit puts to values compatible with zero all parameters except two of 
them, so to reduce the interpolating function to a first order polynomial in /t^. The fit results 
are summarized in Table |2[ We can see that the resulting parameters have a tiny dependence on 
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Figure 4: jScrit vs. p? determined from the susceptibility of chiral condensate, Polyakov loop and plaquette. 
The dashed line is the result of a fit with a polynomial A + Bjx^ of the data at imaginary chemical potential 
(p.^ < 0); the black solid lines enclose the boundary of the uncertainty. The blue vertical line gives the 
position of the first RW transition line. The point in magenta comes from a determination of the peak of the 
susceptibility at j3 = 1.30, taken from Ref. [pl|. 



the observable considered; moreover, the extrapolation of the pseudo-critical line at the jl value 
corresponding to the first Roberge-Weiss transition line, jXnw = iTl/i, is in good agreement with an 



independent determination of the endpoint jSf [|l4|l- This is a confirmation of the structure of the 
phase diagram as sketched in Fig. [I] (right). 

The most important point of the present analysis is to test whether the extrapolation of the 
pseudo-critical line to real jl agrees or not with the two direct determinations of jScrit available so far 
at real jl, i.e. /t=0.20 and 0.30. Such comparison is presented in Figs, ^for each of the observables 
considered. There is an overall substantial agreement; in the case of the chiral condensate there 
might be a deviation at jl=0.30, which calls for a refinement of the numerical analysis. 



4. Conclusions and outlook 



We have presented preliminary results aimed at studying the possibility of the analytic con- 
tinuation of the pseudo-critical line from imaginary to real chemical potential in 2-color QCD. 
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This theory is exempt of the sign problem and, therefore, makes possible to compare analytic con- 
tinuations from imaginary to real chemical potential with direct determinations at real chemical 
potential. 

By determining the position of the peaks in the susceptibilities of three observables (chiral 
condensate, Polyakov loop, plaquette) for varying the temperature at fixed imaginary chemical 
potential, we have interpolated the pseudo-critical Une in the temperature - imaginary chemical 
potential plane. It turns out that the best interpolation for j3crit(M'^) is a first order polynomial in 
jjL^. The extrapolation to real chemical potential of this curve generally agrees with the direct 
determinations at real chemical available so far. A larger statistics and an extension of the data set 
could reveal possible deviations. 
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